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XMM-Newton observation of the persistent Be/NS X-ray 
binary pulsar RX J1 037.5-5647 in a low luminosity state 
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Abstract. The spectra of several X-ray binary pulsars display a clear soft excess, which in most cases can be described with 
a blackbody model, above the main power-law component. While in the high-luminosity sources it is usually characterized by 
low temperature (kT < 0.5 keV) and large emission radius (R > 100 km), in the two persistent and low-luminosity pulsars 
4U 0352+309 and RX J0146.9+6121 this component has a high temperature (kT > 1 keV) and a smaller radius (R < 0.5 
km), consistent with the estimated size of the neutron-star polar cap. Here we report on the timing and spectral analysis of 
RX J1037. 5-5647, another low-luminosity persistent Be binary pulsar, based on the first XMM-Newton observation of this 
source. We have found a best-fit period P = 853.4 ± 0.2 s, that implies an average pulsar spin-up P ~ —2 x 10~ 8 s s _1 in 
the latest decade. The estimated source luminosity is Lx ~ 10 34 erg s -1 , a value comparable to that of the other persistent 
Be binary pulsars and about one order of magnitude lower than in most of the previous measurements. The source spectrum 
can be described with a power law plus blackbody model, with kTsB = 1.26 j^Jg keV and Rbb = 128^21 m > suggesting 
a polar-cap origin of this component. These results strengthen the hypothesis that, in addition to low luminosities and long 
periods, this class of sources is characterized also by common spectral properties. 
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1. Introduction 

Most X-ray binary pulsars (XBPs) are High Mass X-Ray 
Binaries (HMXRBs) in which a Neutron Star (NS) with a mag- 
netic field B ~ 10 12 G accretes matter from a high-mass early- 
type star, either an OB supergiant or a Be star. The X-ray spec- 
tra of both subgroups are generally described by a rather flat 
power law between 0. 1 and 10 keV (photon index ~1) with 
a high-energy cutoff dNagase 2002 ). Several of the brightest 
XBPs (L x ~ 1Q37-38 erg s -i\ fo ave a j S p a m a r ked soft X- 
ray excess dRamsav et ai]|2002l iPaul et alJl2002t iNaik & Paul 



sources with relatively low luminosity (Lx ~ 10 35 36 erg s 
dYokogawa et alj|2000a|: ISasakiet all 120031 : i Haberl & Pietsch 
|2004 [ lUeno et alJ 120041: iMaiid et alJ EooS IHaberl & Pietsch 
l2005h . 



20041 lYokogawa et al]|2000btlKohno et aljboOOtlBurderi et al 
2000t iManousakis et alJ 120091) . A variety of simple models 
(both thermal and non-thermal) have been proposed to fit this 
feature. Only in a few cases this low-energy component unam- 
biguously showed coherent pulses, therefore the debate over its 
origin remains open. 

A large number of new pulsars in BeXRB systems has been 
discovered in the Small Magellanic Cloud, where the observa- 
tions are unaffected by the high interstella r absorption alon 
the Galactic plane dHaberl & Sasaki l2000l: lisrael et all 1200 



Based on the these results. iHickox et alJ d2004l) showed that 
a low-energy component should be visible in all X-ray pul- 
sars that have sufficiently high flux and small absorption. In 
fact, most of the soft-excess sources are at small distances 
and/or away from the Galactic plane (most of them are in the 
Magellanic Clouds). This suggests that the presence of a soft 
spectral component could be a very common, if not an ubiqui- 
tous, feature intrinsic to X-ray pulsars. Different explanations 
for its origin have been proposed, depending on the source lu- 
minosity. When Lx > 10 38 erg s _1 , the luminosity and the 
shape of the soft component can be explained only by the pres- 
ence of an optically-thick accretion disk, which reprocesses the 
hard X-rays coming from the neutron star at its inner edge. In 



ig 
'0; 



less luminous sources, with L 



x 



< 



10 36 erg s" 1 , 



the soft ex- 



Yokogaw a et alj|2003h iMacomb et alj2003l IHaberl et aU2 008 ; 
Eger & Haberll 2008 ). This allowed the detailed study of their 



X-ray spectra down to energies of a few hundred eV, leading 
to the detection of a marked soft excess in most pulsars, also in 



cess can be due to other processes, such as emission by photo- 
ionized or collisionally-heated diffuse gas or thermal emission 
from the surface of the neutron star. Finally, in the sources of 
intermediate luminosity, either or both of these types of emis- 
sion can be present. 
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Recently, a clear thermal exc ess has been observed also 
in the cases of RX J0146.9 +6121 dLa Palombara & Mereghettil 
2006b and 4U 0352+309 dLa Palombara & Mereghettill2007b . 
which are two persistent Be pulsars. This is a specific class of 
binary pulsars characterized by low luminosity (Lx ~ 10 34 ~ 35 
erg s _1 ) and long pulse period (P > 100 s); moreover, they 
show no transient behavi or, since the sourc e luminosity varies 
at most of a factor ~ 10 dNegueruelall 1998b - These properties 
suggest that the Be star orbits the NS in a wide and nearly circu- 
lar orbit, continuously accreting material from the low-density 
outer regions of the circumstellar envelope. This picture has 
been confirmed for 4U 352+309 by the discovery of an or- 
bital period of 250.3 days dDelgado-Marti et al.ll200lb . The de- 
tection of the thermal excess in these two persistent pulsars was 
favoured by their small distance (d < 2.5 kpc) and interstellar 
absorption (7V H ~ 10 21 cm -2 ). 

It is interesting to investigate if other persistent Be binary 
pulsars are characterized by the same type of spectral feature, 
in order to check if it is a common property of this class of 
sources. To this aim, we have observed with XMM-Newton 
the Be/NS binary pulsar RX J 1037. 5 -5647, identified w ith LS 



1698, a B0 Ill-Ve star at - 5 kpc dMotch et all 119971) . This 



sour ce is likely associa ted to the UHURU source 4U 1036- 
56 dFormanet al.ll 1978b . the ARIEL V source 3 A 1036-565 
dWarwick et al. 1 198 lb and the OSO 7 source 1M 1022-554 
(Mar kert et al.1 1979b . Thanks to obser vations performed by 



the PCA instrument on-board RossiXTE, Reig & Roche! (119991) 
performed the first detailed timing and spectral analysis, and 
discovered a pulsation with period P = 860 ± 2 s. The source 
spectrum could be described by a cut-off power law, with a low 
cut-off energy (E cut = 4.7 keV) and a weak (equivalent width 
EQW ~ 65 eV) iron line at 6.5 keV; moreover, the estimated 
source luminosity {Lx — 2 x 10 35 erg s" 1 ) was fully consistent 
with the luminosity level detected by the previous observations, 
and the light curve showed no large variability. Therefore, RX 
J1037.5-5647 is a long-period, low-luminosity Be/NS XBP 
which has been detected every time it has been observed. All 
these pro perties are typ i cal of the persistent Be/NS pulsars, 
therefore iReig & Rochd d 1999b suggested that RX J1037.5- 



5647 is a potential member of this class of sources. 

2. Observations and data reduction 

RX J1037.5-5647 was observed by XMM-Ne wton on 2009 
Janua ry 9. The three EPI C cameras , i.e. o ne pn dStriider et al 



2001b and two MOS dTurner et all 1200 lb . were operated in 
Small Window mode, with a time resolution of 6 ms for the 
pn camera and of 0.3 s for the two MOS cameras; the medium 
thickness filter was used. The effective source exposure time 
was ~ 28 ks for the two MOS and ~ 20 ks for the pn, since 
the use of the Small Window mode for this camera implies a ~ 
30 % dead time. Also the Reflection Grating Spectrometer was 
used but, due to the source faintness, the count statistics was 
very low; therefore we ignored its spectra. 

We used version 8.0 of the XMM-Newton Science Analysis 
System (SAS) to process the event files. After the standard 
pipeline processing, we looked for possible periods of high in- 
strumental background, caused by flares of soft protons with 



energies less than a few hundred keV. We found that the last 
~ 8 ks of the observation were affected by a high soft-proton 
contamination, therefore we rejected the data of this part of 
the observation in the spectral analysis. The resulting effective 
exposure time was 14.5 and 20 ks for the pn and the MOS cam- 
eras, respectively. 

3. Timing analysis 

We selected source counts within a circular region around the 
source position, with extraction radii of 20 and 15" for the pn 
and the two MOS cameras, respectively. Since in Small Window 
mode a very small sky area is imaged by the instrument, we 
were forced to select the corresponding background data on 
small circular areas, with radii of 20". We accumulated all the 
events with pattern range 0-4 (i.e. mono- and bi-pixel events) 
and 0-12 (i.e. from 1 to 4 pixel events) for the pn and the 
two MOS cameras, respectively. This selection resulted in a net 
count rate (CR) of ~ 0.3 cts s _1 for the pn camera and ~ 0. 13 
cts s^ 1 for each of the two MOS cameras. 

To measure the pulse period, we converted the event arrival 
times to the solar system barycenter, for each instrument sep- 
arately; then, in order to increase the count statistic, we com- 
bined the three datasets and mea sured the pulse period by a 
standard phase-fitting technique (IDaU'Osso et alJ l2003b . ob- 
taining a best-fit period P = 853.4 ± 0.2 s. 

Thanks to the XMM-Newton observation, we could inves- 
tigate the pulse profile, for the first time, even at energies be- 
low 3 keV, which could not be studied by RossiXTE. In Fig.[T] 
we show the folded light curves in the three energy intervals, 
0.2-3 (soft), 3-10 (hard), and 0.2-10 keV, together with the 
folded hardness-ratio HR between the hard (H) and soft (S) 
light curves (computed as H/S). The shape of the pulse pro- 
file is similar in the three energy ranges. The measured pulsed 
fraction in the 0.2-10 keV energy band, defined as (CR max - 
CR m i n )/(2 x CR avoragc ), is 75 %. The HR is rather constant 
along most of the pulse period; it increases only at the end of 
the pulse maximum, when the source count rate starts to de- 
crease. 

It is interesting to compare the light curve of the whole 
XMM-Newton observation with the average one, which is ob- 
tained assuming that all the pulses are equal to the average 
folded light curve. To this aim, we considered a binning time 
of 106.5 s, corresponding to eight bins for each pulse period. 
In Fig. |2] we show the total (pn+MOSl+MOS2) background- 
subtracted light curve, in the 0.2-10 keV energy range, and 
the average one, together with their difference. It shows that 
there are significant differences between the individual pulses, 
as well as some variability on longer timescales. Moreover, we 
checked that the source HR between the hard and soft bands 
is rather constant around the average value of 1.08 (x 2 /d.o.f. 
= 1.89/134), and does not depend on the source count rate 
(Fig.E). 

4. Spectral analysis 

For the source and background spectra, we considered the same 
extraction parameters we used for the light curves. We gener- 
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Fig. 2. Top: observed (black) and average {red) light curve of RX J1037. 5-5647 in the energy range 0.2-10 keV; Bottom: count- 
rate difference between the individual time bins of the observed and average light curve. 
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Fig. 1. Background-subtracted light curves of RX J1037.5- 
5647 in the energy ranges 0.2-3, 3-10, and 0.2-10 keV, folded 
at the best-fit period P = 853.4 s. The vertical lines indicate the 
phase intervals used for the phase-resolved spectroscopy (see 
Section 5). 



ated the applicable response matrices and ancillary files using 
the SAS tasks rmf gen and arf gen. In order to ensure the 
applicability of the \ 2 statistics all spectra were rebinned with 
a minimum of 30 counts per bin; they were fitted in the en- 
ergy range 0.2-10 keV using XSPEC 1 1.3.2. In the following, 
all spectral uncertainties and upper-limits are given at 90 % 
confidence level for one interesting parameter. 

After checking that separate fits of the three cameras gave 
consistent results, we fitted them simultaneously, in order to in- 
crease the statistics; to this aim, we introduced relative normal- 
ization factors among the spectra of the three cameras. Using 
an absorbed power-law (PL) model, we obtained a hydrogen 
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Fig. 3. Hardness-ratio variation of RX J 1037.5-5647 as a func- 
tion of the 0.2-10 keV count rate, with a time bin of 213 s. 



column density Nn = (6.2 ± 0.5) x 10 21 cm 2 and a photon- 
index T = 0.94±0.03, with x^/d.o.f. = 1.30/277. On the other 
hand, using an absorbed blackbody (BB) model, we obtained a 
hydrogen column density N B = (0.9 ± 0.2) x 10 21 cm~ 2 and 
a BB temperature kT BB = 1.52 ± 0.03 keV, with x^/d.o.f. = 
1.20/277; assuming a source distance of 5 kpc we obtained a 
BB radius R B b = 131 ± 4 m. 

In both cases the fit of the spectrum with a single- 
component model shows large residuals, therefore we repeated 
the fit with a PL+BB model. In this way we obtained a sig- 
nificant improvement of the fit quality (Fig. |4j, since we ob- 
tained x^/d.o.f = 0.88/275; the corresponding best-fit param- 



eters are iV H = (2.8 ± 0.9) 



T = 0.51 



+0.17 
-0.29 



and A;Tbb = 1.26_q'qq keV. In comparison with the single PL 
and BB model, the F-test analysis provided a probability, re- 
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spectively, P = 1.26 X 10" 24 and P = 9.95 x 10" 20 that the 
improvement of the fit occurs by chance. Both components are 
significant at 99 % confidence level: the PL normalization is 
ip L = 7.l±fl x 1(T 5 ph cm" 2 s- 1 keV" 1 at 1 keV and, 
assuming a source distance of 5 kpc, we obtained a radius 
i?BB = 128^21 m f° r tne BB component. The unabsorbed flux 
in the energy range 0.2-10 keV is fx ~ 4 x 10 



12 — ? 

erg cm 



s , about 42 % of which is due to the BB component; this 
translates into a source luminosity Lx — 1.2 x 10 34 erg s _1 . 
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Fig. 4. Top panel: total spectrum of RX J1037.5-5647 with the 
best-fit PL+BB model. The spectra of the pn, MOS1 and MOS2 
cameras are shown in black, red and green, respectively. Middle 
panels: residuals (in units of er) between data and model in the 
case of the single PL and of the single BB. Bottom panel: resid- 
uals in the case of the PL+BB. 

We looked also for narrow iron K a emission lines between 
6 and 7 keV, with different widths between and 0.5 keV. We 
found no evidence for such a component, with an upper limit 
on its equivalent width of 0.2 keV (at 90 % c.l.) at most. 

5. Phase-resolved spectroscopy 

In the source folded light curve (Fig. Q]) we observed that the 
source spectrum becomes slightly harder at the end of the pulse 
maximum. In order to study the source behavior in more detail, 
we analyzed the background subtracted spectra in two different 
phase intervals, i.e. <f> = 0-0.68, 0.92-1 (phase interval A) and 
4> = 0.68-0.92 (phase interval B) 

We fit the two background-subtracted spectra with the 
best-fit PL+BB model of the phase-averaged spectrum, leav- 



ing only the relative normalization factors between the two 
phases free to vary; in Fig. [5] we report the ratios of the two 
spectra of each instrument to this renormalized average model. 
The data of the three detectors show that the spectra of phase 
intervals A and B are, respectively, softer and harder than the 
phase-averaged one, which implies a spectral variability along 
the pulse phase. 




Energy (keV) 

Fig. 5. Ratios of the spectra corresponding to the phase inter- 
vals shown in Fig. Q] to the average model spectrum. The pa- 
rameters are fixed at the best-fit values for the total spectrum, 
except for the normalizations, which are at their best-fit values 
for each phase interval. 

We then fitted the two spectra independently, using both a 
single PL model, a single BB model and a PL+BB model; their 
best-fit parameters are reported in Table [T] In both cases we 
found that, as for the phase-averaged spectrum, the BB model 
provides a better fit than the PL model and that the use of a 
PL+BB model improves significantly the spectral fit goodness, 
since there is a large decrease in the \ 2 value and the F-test 
analysis gives a negligible probability that the improvement of 
the fit occurs by chance. Irrespective of the adopted model, we 
always found that the absorbing hydrogen column density is 
consistent with a constant value, as expected if it is due to the 
interstellar extinction. With both the single-component and the 
two-component models the photon index shows evidence for a 
variability between the two phase intervals, thus confirming the 
spectral variability; on the other hand, the BB temperature and 
radius are variable only in the case of the single-component 
model, while they are rather constant in the case of the two- 
component model, since there is a full agreement for the emis- 
sion radii and also the temperatures are comparable. 

To investigate the relative variations of the two components 
with the period phase, we also simultaneously fitted the two 
spectra forcing the same value of Ah, T and kT-QB for the two 
phase intervals. In this case we obtained A H = l-4±gj x 10 21 
cm" 2 , r PL = -0.2t°;3, and kT BB = 1.21±g^ keV, with 
X 2 /d.o.f. = 0.99/115; the corresponding normalization values 
are reported in Table [2] With these constraints on the model 
parameters, the spectral changes as a function of the phase are 
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Table 1. Best-fit spectral parameters for the phase-resolved spectroscopy of RX J1037. 5-5647, in the case of the independent 
fit of the two spectra. Nn, /cTbb, and Rbb are measured in units of 10 21 cm -2 , keV, and meters (for a source distance of 5 kpc), 
respectively. Pp_ test (PL) and Pp—test [BB) are the probability that the improvement of the fit, compared to the single PL and 
BB model, respectively, occurs by chance. 



Phase interval 


A 


B 


A 


B 


A 


B 




6.6 ± 0.7 


6.6±i;g 


0.5 ± 0.3 


0+u.s 

"-0.0 


2 l +i ' 4 


O o+u.y 


r 


1.01 ± 0.06 


0.79 ±0.10 






Q4+U.6 
u/f -1.0 


9 c+0.3 

" z - J -0.2 


JPL 


4.8 (100%) 


5.5 (100 %) 






2.0 (51 %) 


1.5(31 %) 


Rbb 






1.45 ± 0.03 
139 ±6 


1171? 


i 99+U.I8 
L - z - z --0.09 

142+29 


i.44l^;- 
uitll 


/bb 






3.4 (100 %) 


4.1 (100%) 


1.9 (49 %) 


3.4 (69 %) 


/tot 


4.8 


5.5 


3.4 


4.1 


3.9 


4.9 


*2/«Lo.f. 


2.00/80 


1.50/33 


1.45/80 


1.17/33 


1.05/78 


0.80/31 


P F -test(PL) 










4.79xl0~ 12 


2.15xl0" 5 


P F -test(BB) 










1.33xl0" 6 


LlOxlO -3 



a Unabsorbed flux in the energy range 0.2-10 keV, in units of 10 1 erg cm s 



reproduced by the variations in the relative contribution of the 
two components. The values reported in Table |2] show that al- 
most all of the spectral variability is due to the PL component, 
since its flux changes of a factor ~ 2 between the two phase 
intervals; on the other hand, there is only a 15 % variation of 
the BB flux, since we obtain comparable values for the emis- 
sion radius. However, this finding is affected by large errors on 
the normalization of the PL component. 

As further check, we modified the test model by linking to- 
gether the BB parameters for the two spectra, while both the 
photon-index V and the normalization Zpl of the PL com- 
ponent could vary independently in the two phase intervals. 
The resulting best fit gives /V H = (2.3 ± 1.0) X 10 21 cm" 2 , 
fcT BB = 1.25±g;£jj keV and R BB = 137±H m, while the 
power-law parameters are shown in Table |2 Even with this 
model we found a very good fit quality (x 2 /d.o.f. = 0.96/1 15), 
fully comparable to that obtained in the previous case. This 
result suggests that it is possible to attribute the whole spec- 
tral variability to the PL component. However, we note that the 
power-law parameters are not well constrained, since they are 
always affected by large errors. Moreover, for the spectrum of 
the phase interval B the independent fit with a PL+BB model 
provides a large negative value of T, which is very different to 
that obtained from any other fit. This result can be explained 
by the large contribution of the BB component to the total flux, 
while in the other cases most of the flux is due to the PL com- 
ponent: in fact, for the phase interval B we can obtain a good 
fit also using a simple BB model (x 2 /d.o.f. = 1.17/33). This 
means that an additional PL component is required only to fit 
the high energy end of the spectrum, thus providing a large 
negative value of the photon index. 



6. Discussion 

Before the XMM-Newton observation of RX J1037.5-5647 
the only X -ray investigation of this source was performed by 
RossiXTE (Reig & R ochelll999h . Therefore it is interesting to 
compare our results with those obtained ten years before. 



Table 2. Best-fit values for the black-body and power-law pa- 
rameters, when the two spectra are fitted simultaneously with 
common values of 7Vh, T and kT B B (case 1) or with common 
values of /Yh, kT B B and Rbb (case 2). 



Phase 
interval 


case 1 
A B 


case 2 
A B 






2.3 ±1.0 2.3 ±1.0 


r 

-'PL 
fpL 


n 9+U.3 n 9+U.3 
""' z -0.9 " lAZ -0.9 
, 9+2.0 9 4+3.2 
i - z -1.0 z "*-2.0 

1.5(39%) 2.9(60%) 


05 ,+u.3b 00 o+'J.37 

u - J1 -0.57 u - u ^-0.66 

60 +5.i 43+ 41 

O.U_ 4 g ^• J -3.2 

2.0(51%) 2.9(61%) 


fcT BB (keV) 

%B 

/bb 


, 9,+U.UV , 9,+U.UV 
J - z -0.06 J - ZJ -0.06 

mtl 1 i47l^ 

2.3 (61 %) 2.0 (40 %) 


; 9<r+U.ia , 9<-+U.13 
-'• z -'-0.09 '■ ZJ -0.O9 

137±H 137tfx 
1.9(49%) 1.9(39%) 


/tot 
X 2 /d.o.f. 


3.8 4.9 
0.99/115 


3.9 4.8 
0.96/115 



a Hydrogen column density in units of 10 21 cm 2 

b Intensity of the power-law component in units of 10 -5 ph cm" 2 s" 1 

keV" 1 at 1 keV 

c Unabsorbed flux in the energy range 0.2-10 keV, in units of 10 -12 
erg cm" 2 s -1 

d Radius of the black-body component (in metres) for a source dis- 
tance of 5 kpc. 

From the timing analysis, we have obtained a new refined 
pulse period, P = 853.4 ± 0.2 s, which is lower than the 
previous value of 860 ± 2 s found by RossiXTE. This im- 
plies a pulsar spin-up during the past ten years, with P = 
(—1.91 ± 0.58) x 10" 8 s s _1 , and suggests a momentum trans- 
fer to the neutron star. 

Based on the observed source flux, we have estimated a 
source luminosity Lx ~ 10 34 erg s _1 in the 2-10 keV en- 
ergy range (assuming a source distance of 5 kpc), which is one 
order of magnitude lower than the luminosity level detected 
by RossiXTE. From the spectral analysis, we have obtained 
a hydrogen column density A^h = (2.8 ± 0.9)xl0 21 cm" 2 , 
which is much lower than the values estimated by RossiXTE: 
the lowest result (obtained with a double BB model) was (2.1 
± 0.3) xlO 22 cm" 2 , which is almost one order of magni- 
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tude higher than our result; moreover, using our same PL+BB 
model, the RossiXTE data provided a best-fit value of (8.2 ± 
0.5)xl0 22 cm" 2 , which is ~ 30 times higher than our value. 
On this subject, we emphasize that the energy range of the 
RossiXTE spectral analysis (above 3 keV) is not well suited for 
a good estimate of /Vh, while the low-energy end of the XMM- 
Newton spectra (0.2 keV) allows a much more reliable analy- 
sis. Moreover, we also note that LS 1698, the optical coun- 
terpart of RX J103 7.5-5647, has a color excess E(B-V) ~ 0.75 
(Motch e t al.l 1997 ): assuming Ay = 3.1 E(B-V) and the average 
relation Ay = Nu x 5.59 x 10 ~ 22 cm" 2 between optic al ex- 
tinction and X-ray absorption (Pred ehl & Schmittlll995l) . this 



would predict TVjj = 4.16 x 10 21 cm" 2 , a value comparable 
to our result. This indicates that in the XMM-Newton observa- 
tion the measured absorption is due to the interstellar medium 
and not to local matter around the system, which has a negli- 
gible effect, while in the RossiXTE observations also intrinsic 
absorption is possible. 

The RossiXTE spectrum analyzed by Reig & Roche (1999) 
was fitted with the same model {PL+BB), resulting in a quite 
steeper power-law component (T = 1.84 ± 0.06 instead of 
0.51±°;^) and in a blackbody component with a much higher 



1 26+ 016 



temperature (kT = 2.9 ± 0.2 keV instead of kT 
keV). Even in this case, it is possible that these differences arise 
not only because the spectral analysis was performed in two 
different energy ranges, but also because the source was one 
order of magnitude brighter during the RossiXTE observation. 
On the other hand, the two analysis provided comparable val- 
ues for the emission radius of the blackbody component (100 
± 10 m for RossiXTE and 128± 2? m for XMM-Newton). 

We found no evidence of an Fe-K line in the source spec- 
trum, with an up per limit of 0.2 keVto its equivalent width (at 
90 % c.l.), while lReig & Rochd (fl999) report a gaussian line at 
E = 6.5 ± 0.2 keV, with an equivalent width ~ 65 eV. It is very 
likely that XMM-Newton missed a clear detection of this fea- 
ture due to the different luminosity level of the source, since a 
lower accretion rate on the NS implies a smaller amount of ma- 
terial around it and, then, a reduced probability for the emerg- 
ing X-rays to produce an iron line feature. 

The comparison of the pulse profiles in two energy bands 
suggests a possible spectral variability of the source with the 
pulse phase, therefore we performed a phase-resolved spectral 
analysis, which previously had never been carried out for this 
source. The fit of two phase-resolved spectra with the best- 
fit model of the phase-averaged one confirmed the suspected 
variability. The independent fit of these spectra with different 
emission models showed that the absorbing hydrogen column 
density is nearly constant along the pulse phase, as expected 
if it is due to the interstellar extinction; therefore the observed 
variability must be due to the intrinsic source emission. The 
same fit proved also that, for both spectra, both the thermal and 
non-thermal component are necessary in order to well repro- 
duce the observed data. Finally, the simultaneous fit of the two 
spectra with a common value of iVn shows that the whole spec- 
tral variability can be attributed to the PL component, since the 
BB one is nearly constant along the pulse. 



Reig & Rochel d 19991) suggested that RX J1037.5-5647 is 
a persistent Be binary pulsars, since this source has been de- 
tected every time it has been observed and does not show large 
outbursts. In order to investigate the long-term evolution of RX 
J1037. 5-5647, we considered all the X-ray flux measurements 
obtained for this source since its disc overy in t he early '70s. 
For th e UHURU dForman et alJll978l). O SO 7 dMarkert et al ' 



1979) and ARIEL V dWarwick et al.lll98ll) detections, we used 



the published CR values and the appropriate count-rate-to- 
flux conversion factors to obtain the corresponding 2-10 keV 
fluxes. In the case of the RossiXTE All Sky Monitor (ASM), 
we considered the day-by-day light curve, that covers the time 
interval from January 1996 to March 2009. Since the detec- 
tion sensitivity of this instrument is very low, in several cases it 
failed to detect RX J1037. 5-5647, owing to its low flux: for this 
reason, in order to estimate the average flux level of the source, 
we binned the data in five different time bins of 1000 days 
each; we calculated the average CR in each bin and, assuming 
our best-fit PL+BB emission model, we used the HEASARC 
tool WebPIMMS to obtain the corresponding average flux be- 
tween 2 and 10 keV. On the other hand, for the data of the 
RossiXTE Photo n Counting Array (PC A), we based on the re- 
sults obtained bv lReig & Roche ( 1999) to calculate the source 
flux in the same energy range. Finally, we also considered the 
CRs measured by ROSAT, which were obtained during the All 



Sky Survey ( Voges et a 



servation ( Mot ch et al 



19991) and in a following, specific ob- 



1997). In their case, we had to infer the 



source flux in the 2-10 keV energy band from the CR mea- 
sured in the 0.1-2.4 keV range: to this aim, we used again 
the WebPIMMS tool and assumed our best-fit PL+BB emission 
model. 



Assuming a source distance of 5 kpc (Mo tch et alJ l!997). 
we used the calculated fluxes to estimate the corresponding lu- 
minosities, which are reported in Fig. [6] It shows that, in most 
observations performed before XMM-Newton, RX J1037.5- 
5647 was detected at a luminosity level Lx = (1—3) x 10 35 erg 
s" 1 between 2 and 10 keV; moreover, also the average luminos- 
ity observed by the ASM is ~ 10 35 erg s" 1 . This is a remark- 
able behavior, since the performed observations span a time in- 
terval of about 35 years, and confirms that RX J1037. 5-5647 is 
a persistent, low-luminosity binary pulsar. However, the obser- 
vation performed by XMM-Newton detected the source at a lu- 
minosity level Lx ~ 10 34 erg s _1 , which is one order of mag- 
nitude lower than the previous average value. Moreover, in one 
of the ROSAT observations it was detected with an even lower 
luminosity. It is possible that the ROSAT results are due to an 
inaccurate estimate of the source flux, since it is obtained from 
a CR measurement in a different energy band, but even such 
an error could hardly explain the very low luminosity value ob- 
tained in the second ROSAT observation. These results confirm 
that RX J1037. 5-5647 is a persistent source, since it has been 
detected every time it has been observed, but show also that the 
source is characterized by a significant variability. 

A BB excess above the main power-law compo- 
nent has been detected also i n the case of LS I +61° 
235/RX J0146.9+6121 dLa Palombara & Mereghetti 



2001 



| 2006l) and X Persei/4U 0352+ 309 dCoburn et al 
lLa Palombara & Mereghettil 20071) . which are the two main 
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Fig. 6. Luminosity history of RX J1037.5-5647 from 1971 to 
2009. Reported luminosities are based on the una bsorbed flux 



in the 2-10 keV and on a source distance of 5 kpc dMotch etal 



1997b . Open circles refer to the average luminosity (between 2 
and 10 keV) estimated by RossiXTE ASM in five time intervals 
of 1000 days each. 



members of the class of the persistent Be/NS XBPs. A 
similar feature was observed also during the low-luminousity 
quiescent state of the transient Be/NS XBP 3A 0535+262 



dOrlandini et all l2004t iMukheriee & Paull 12005b . which is 



comparable to the previous persistent Be pulsars for its low 
luminosity (Lx < 10 35 erg s _1 ) and long pulse period (P = 
103 s); moreover, a blackbody component has been observed 
also in 4U 2206 +54 dMasetti et al.ll2004t ITorrejon et alll2004 ; 
Reig etal 2009), a high-mass X-ray binary with a comparable 
luminosity and a very long pulse period (P = 5560 s), and in 
the X -ray pulsar (P = 358.6 s) SAX J2103.5+4545 dinam et al 



2004). In Fig. [7] we report the best-fit radius and temperature 
for these sources, together with lines showing four different 
levels of the blackbody luminosity. In the case of RX J1037.5- 
5647 (filled circles) and 4U 0352+309 (filled squares), we 
report two set of values, since both sources were observed by 
RossiXTE and XMM-Newton at two different luminosity levels 
(Lx ~ 10 34 and ~ 10 35 erg s _1 , respectively): in the case of 
RX J1037. 5-5647 the low and high luminosity were observed 
by XMM-Newton and RossiXTE respectively, while for 4U 
0352+309 it was the opposite. We report various measurements 
also for 4U 2206+54 (crosses), corresponding to different 
observations. The figure shows that the results obtained for 
RX J1037. 5-5647 are in full agreement with those obtained 
by various observations of the previous sources, since their 
spectral parameters are clustered in a narrow range of values, 
i.e. /cTbb ~ 1-1-5 keV and Rbb ~ 100 m. We emphasize 
that, in all these cases, the estimated total source X-ray 
luminosity is ~ 10 34 erg s _1 , with a 20-40 % contribution 
of the blackbody component. For the persistent Be pulsars 
RX J1037.5-5647 and 4U 0352+309, spectral values outside 
the previous range (points 2 and 4) were obtained by the 
XMM-Newton observation of 4U 0352+309 and the RossiXTE 
observation of RX J1037. 5-5647, when the sources were 
detected at a higher luminosity level (Lx ~ 10 35 erg s" 1 ); 



however, we note that comparable blackbody parameters were 
obtained also for SAX J2103. 5+4545, which is characterized 
by a comparable luminosity. On the other hand, in most cases 
the non-Be binary pulsar 4U 2206+54 shows a different 
behavior also when observed at a comparable luminosity level: 
kT B B < 1 keV and Rbb > 500 m (points 8 to 11). Finally, 
we note that the XMM-Newton spectrum of RX J 1037. 5-5647 
is characterized by the hardest PL component among the 
previous sources, since it has a photon-index T = 0.5, while 
r > 1 in all the other sources. 
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Fig. 7. Best-fit values for radius and temperature of the 
BB component in the case of RX J1037. 5-5647 (filled cir- 
cles), 4U 0352+309 (filled squares), RX J0146.9+6121 (open 
circle), 3A 0535+262 (open square), SAX J2103.5+4545 
(open triangle), and 4U 2206+54 (crosses). The contin- 
uous lines connect the blackbody parameters correspond- 
ing to four different levels of luminosity (i n erg s _1 ). 
References: 1 - this work; 2 - |R eig & Roche (199 9J); 3 - 



Coburnetal. (2001); 4 - lLa Palom bara & Mereghetti l d200 



5 - IMukheriee & Paull d2005l): 6 - l La Palombara & Mereghett; 



d2006h: 7 -iMasetti et alJd2004t): 8.9.10 -ITorrejon et alJ d2004: 



11 - iReig et alJ d2009h : 12- linam etalld2004 . 



An X- ray excess has been observed also in several other 
XBPs (see lLa Palombara & Mereghettill2006l for a review) but, 
contrary to the previous sources, in their case the fit of this ex- 
cess with a thermal emission model provided low temperatures 
(kT < 0.5 keV) and large emitting regions (R > 100 km); for 
this reason, this feature is usually described as a soft excess. In 
Fig. [8] we report the luminosity and pulse period of the XBPs 
with a detected thermal excess. They are divided in two well 
distinct groups: the sources in the first group are characterized 
by high luminosity (Lx > 10 37 erg s _1 ) and short pulse pe- 
riod (P < 100 s), and in most cases they are in close binary 
systems with an accretion disk; those in the second group have 
low luminosities (Lx < 10 36 erg s _1 ) and long pulse periods 
(P > 100 s), since they have wide orbits and are wind fed sys- 
tems. Among the sources of the second group of XBPs, the six 
pulsars discussed above (reported as asterisks and filled circle) 
are the ones that have, at the same time, the lowest luminosi- 
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ties and the longest periods with the only exception of SAX 
J2103. 5+4545, which has the highest luminosity of this group 
of sources. They are characterized by a BB component with a 
high temperature and a small emission radius, while the others 
show a soft thermal excess with low temperature (kT < 0.5 
keV) and large emission area (R ~ a few hundred km). The 
hot BB spectral component separates these low-luminosity and 
long-period sources from all the other pulsars, strongly sug- 
gesting that they form a distinct and well defined class of binary 
pulsars. 



□ □ 



0.01 0.1 1 10 100 1000 10 4 

Pulse Period (s) 

Fig. 8. X-ray luminosity (in the 2-10 keV energy range) of the 
XBPs with a detected thermal excess as a function of the pulse 
period. The filled circle refers to the XMM-Newton observation 
of RX J1037. 5-5647, asterisks refer to other detections of hot 
BB pulsars, open squares to all the other bright sources. 

Based on ou r resul ts and on the emission models proposed 
by lHickox etail d2004 . the luminosity of RX J1037.5-5647 is 
too small for an interpretation of the soft excess in terms of re- 
processing of the hard X-ray photons in optically thick accret- 
ing material. Moreover, both the thermal emission and the hard 
X-ray reprocessing in a diffuse, optically thin gas around the 
neutron star are unlikely, since these processes would not give 
a blackbody spectrum. We therefore favor the interpretation of 
the soft excess in RX J 1037. 5-5647 as thermal emission from 
the neutron star polar cap. If we assume that the source is in the 
'accretor' status, with matter accretion on the NS surface, the 
blackbody emitting radius of ~ 130 m is consistent with the ex- 
pected size of the polar cap. In fact, assuming Mns = 1.4 M 
and i?NS = 10 6 cm, the source luminosity Lx = 1.2 x 10 34 erg 
s _1 implies an accretion rate M ~ 6.4 x 10 13 g s _1 and, adopt- 
ing i?NS = 10 12 G, a mag netospheric radius R m ~ 2.5 x 10 9 
cm ( Campana et al.|[l998l) . In this case, based on the relation 
Rcoi ~ i?NS (i?Ns/#m) ' 5 faickox et alJ2004h . we would ob- 
tain R co i ~ 200 m, comparable to the estimated blackbody 
emitting radius 

The above results suggest that the observed thermal com- 
ponent can be attributed to the emission from the NS polar 
caps, as already proposed for 4U 0352+309, RX J0146.9+6121 
and 3A 0535+262. If this description is correct, we would 
expect to observe some variability of the thermal component 



along the pulse phase. This hypothesis cannot be confirmed by 
the present data, since the phase-resolved spectroscopy shows 
that the significant spectral variability detected along the pulse 
phase can be due to the PL component, with no variation of the 
BB one. 



7. Conclusions 

We have reported the analysis of the data collected by XMM- 
Newton in a ~ 28 ks observation of the Be/NS X-ray pulsar RX 
J1037. 5-5647. The unabsorbed flux corresponds to a source 
luminosity Lx — 1.2 x 10 34 erg s _1 in the 2-10 keV energy 
range, which is about one order of magnitude lower than the 
average level of the previous observations. The source pulse 
period is P = 853.4 ± 0.2 s, which, together with previous 
measurements, implies an average pulsar spin-up P ~ — 2 x 
10 -8 s s _1 in the latest decade. 

We could perform the first accurate timing and spectral 
analysis below 2 keV for this source. The pulse profile shows 
a complex structure, which is not sinusoidal but does not de- 
pend on energy. In the source spectrum we detected a count 
excess above the main power-law component, well described 
by a blackbody with kT-QB — 1.3 keV and Rbb — 130 m, 
which contributes for ~ 40 % to the total source luminosity. 
We found no evidence of a narrow iron K a line between 6 and 
7 keV, with an upper limit of 0.2 keV on its equivalent width. 

The new XMM-Newton detection of the long-period Be bi- 
nary pulsar RX J1037. 5-5647 at a low luminosity level con- 
firms that it belongs to the class of the persistent Be binary 
pulsars, whose main representatives are 4U 0352+309 and RX 
J0146. 9+6121. Moreover, there is a full agreement among the 
results of the spectral analysis of these three sources; there- 
fore it is reinforced the hypothesis that a thermal excess of 
blackbody type, with high temperature {kT > 1 keV) and 
small emission area (R < 0.5 km), is a common feature of this 
type of sources. However, this type of excess is not exclusive 
of the persistent low-luminosity Be binary pulsars, since also 
the transient source 3A 0535+262 shows this component. This 
means that, on one hand, all the low-luminosity Be X-ray bi- 
nary pulsars show this component but, on the other hand, the 
presence of such component in a Be binary does not imply that 
the system is a persistent low-luminosity source. 

As in the other two sources, also for RX J1037. 5-5647 the 
blackbody radius is comparable to the estimated size of the NS 
polar-cap, which suggest that the origin of this component is 
on the NS surface. However, the phase-resolved spectroscopy 
does not confirm this scenario, since the observed spectral vari- 
ability along the pulse period can be attributed to the power- 
law component and is consistent with a constant thermal com- 
ponent. 
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